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Abstract 

A study of the effects of meson-exchange currents and isobar configurations in A{e, e'p)B 
reactions is presented. We use a distorted wave impulse approximation (DWIA) model where 
final-state interactions are treated through a phenomenological optical potential. The model 
includes relativistic corrections in the kinematics and in the electromagnetic one- and two-body 
currents. The full set of polarized response functions is analyzed, as well as the transferred 
polarization asymmetry. Results are presented for proton knock-out from closed-shell nuclei, for 
moderate to high momentum transfer. 
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1 Introduction 



For the last decades coincidence (e,e'p) reactions on complex nuclei have provided precise 
information on bound nucleon properties, which have made it possible to test carefully 
the validity of present nuclear models [1, 2, 3, 4]. Although the analysis of these processes, 
making use of different distorted wave approaches and coupled-channel models, has been 
extremely useful, there are still uncertainties associated to the various ingredients that 
enter in the description of the reaction mechanism: treatment of final-state interactions 
(FSI), nuclear correlations, off-shell effects, Coulomb distortion of the electrons, relativis- 
tic degrees of freedom, meson-exchange currents (MEC), etc. All of these ingredients 
affect the evaluation of the differential cross section and hence lead to ambiguities in the 
extraction of the spectroscopic factors. The origin of this uncertainty is directly connected 
with the complexity of the dynamics of the reaction and the different approaches to han- 
dle it, which produce different cross sections. It is clear that a reliable determination of 
spectrospic factors requires an accurate description of the reaction mechanism. Important 
efforts in this direction have been done in recent works [5, 6, 7, 8, 9]. 

The measurement of the separate nuclear response functions and asymmetries imposes 
additional restrictions over the theory. The exclusive response functions, which include 
different components of the hadronic tensor taken along the longitudinal (L) or transverse 
(T) directions with respect to the momentum transfer q, may present very different sensi- 
tivities to the different aspects of the reaction. In this sense, it is interesting to point out 
that MEC are shown to contribute mainly to the transverse components [10, 11, 12], while 
relativistic degrees of freedom play a crucial role in the interference TL response [8, 13]. 
Thus, a joint analysis of cross sections and response functions, comparing the experimental 
data with the theoretical predictions, can provide very relevant and complementary infor- 
mation on the reaction mechanism. Separate response functions and the TL asymmetry 
have been measured for 16 0(e,e'p) at moderate [14, 15] and high [16] g-values. The asym- 
metry A TL , obtained from the difference of cross sections measured at opposite azimuthal 
angles (with respect to q) divided by the sum, results particularly relevant because it 
does not depend on the spectroscopic factors. For high missing momentum p > 300 
MeV/c, Atl presents an oscillatory structure that has been shown to be consistent with 
predictions of 'dynamical' relativistic calculations [13, 17, 18, 19]. 

The advent of longitudinally polarized beams [20] and recoil polarization measure- 
ments [21, 22] has importantly enlarged the number of observables which can be accessi- 
ble with this type of experiments, a fact that is welcome to challenge the theory strongly. 
In recent experiments carried out at MIT-Bates and Jefferson Lab, the induced (P) and 
transferred (P') polarization asymmetries were measured for complex nuclei, 12 C [21] and 
16 [22], respectively. In both cases (q, c<j)-constant kinematics has been selected with 
q m 760 MeV/c, u m 290 MeV at MIT-Bates and q m 1000 MeV/c, u « 450 MeV 
at TJlab. Since the transfer momentum values are high enough, relativistic degrees of 
freedom should be incorporated in a consistent description of these reactions. After the 
pioneering work in [23, 24], a detailed study on the induced normal polarization P n has 
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been presented in [25, 26] within the framework of the relativistic distorted wave impulse 
approximation (RDWIA). A comparison with non-relativistic analyses was also discussed. 
The sensitivity of polarized observables to channel coupling in final state interactions was 
analyzed in [3, 27], while in [17] the study was focussed on the effects of spinor distortion 
over the transfer polarization ratio P' x /P' z - In [28] the whole eighteen recoil nucleon polar- 
ized responses were computed from intermediate to high momentum transfer in the Dirac 
eikonal formalism. A comparison between the predictions of the Glauber and eikonal 
models for P n was presented in [29] with the aim of bridge the gap between the low- and 
high-energy description of FSI. More recently a theoretical study of kinematical and dy- 
namical relativistic effects over polarized response functions and polarization asymmetries 
has been performed in [30, 31] within the relativistic plane-wave impulse approximation 
(RPWIA). A general analysis of all the polarized observables within the RDWIA is at 
present in progress [32]. 

Our main aim in this work is to explore in depth the role played by the two-body 
currents in recoil nucleon polarization observables. Some previous analyses on this subject 
have been done by the Pavia group [33, 34] and the Gent group [35, 36]. The calculation 
of MEC in [33, 34] makes use of an effective one-body operator leading to results which, in 
the unpolarized case, differ significantly from those obtained with other approaches that 
describe properly the two-body currents [37, 38]. Recently the unpolarized model of [34] 
has been improved in [39], but differences with other MEC calculations still persist [40]. 
In [36] the induced and transferred polarization asymmetries P n , P{ and P[ were evaluated 
for different kinematical situations. The model considered did not rely on any empirical 
input with respect to the FSI, describing the bound and scattering states as the solutions 
of the Schrodinger equation with a mean field potential obtained from a Hartree-Fock 
calculation. MEC were included based on the formalism developed in [37] which also 
differs from the MEC analysis performed in [38, 40]. In addition, in [36] results for high 
momentum transfer (up to q — 1 GeV/c) were evaluated including relativistic corrections 
into the one-body current operator obtained through the Foldy-Wouthuysen method. 

In this work we extend the DWIA+MEC model developed for unpolarized reactions 
in Refs. [38, 40] in order to describe the spin observables in (e, e'p) processes from closed 
shell nuclei. This model takes care of relativistic degrees of freedom by making use of 
semi-relativistic (SR) operators for the one-body (OB) current [13, 41, 42, 43] as well 
as for the two-body MEC [40, 44, 45, 46]. The SR currents are obtained by a direct 
Pauli reduction of the corresponding relativistic operators by expanding only in missing 
momentum over the nucleon mass while treating the transferred energy and momentum 
exactly. Relativistic kinematics for the ejected nucleon is assumed throughout this work. 
Finally FSI are incorporated through a phenomenological optical potential which, for high 
momentum transfer, is taken as the Schrodinger-equivalent form of a S-V Dirac optical 
potential. The goal of this work is to use the SR approach to evaluate the importance of 
MEC effects upon the spin observables and their dependence on the FSI for intermediate 
to high momentum transfer. As a complete relativistic distorted wave analysis of MEC 
in (e,e'p) processes is still lacking — the only study in this direction has been performed 
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taking into account only the contact current [47] — the use of the SR model becomes, 
as a starting point, a convenient way of implementing relativistic effects in existing non 
relativistic descriptions of the reaction mechanism in order to explore the high momentum 
region. 

The paper is organized as follows: in Section 2 we outline the DWIA formalism de- 
scribing in detail the multipole expansion done for the separate response functions. In 
Section 3 we present our results for the polarized response functions and transferred po- 
larization asymmetries for selected kinematics near the quasielastic peak. Finally our 
conclusions are drawn in Secion 4. 



2.1 Cross section and response functions 

The general formalism for coincidence electron scattering on nuclei involving polarization 
degrees of freedom has been presented in detail in Refs. [4, 23, 24, 48]. In this section we 
simply provide the basic description of our DWIA model focusing on the development of 
the multipole expansion used to compute the response functions. For this end we follow 
closely the multipole formalism developed in [49] for polarized nuclei. 

We consider the process depicted in Fig. 1, in which an incident electron with four- 
momentum K£ = (e e , k e ) and helicity h interacts with a nucleus A, scatters through an 
angle 6 e to four-momentum K'^ = (e', k' e ) and is detected in coincidence with a nucleon 
with momentum p' and energy E'. The four-momentum transferred to the nucleus is 

= K^ — K 1 ^ = (ou, q), verifying Q 2 = u 2 — q 2 < 0. The polarization of the final nucleon 
is measured along an arbitrary direction defined by the unitary vector s. Assuming plane 
waves for the electrons and neglecting the nuclear recoil, the cross section can be written 
in the extreme relativistic limit (ERL) m e < e e , as [48] 



where a separation has been made into terms involving polarized and unpolarized incident 
electrons. Using the general properties of the leptonic tensor it can be shown that both 
terms, £ and A, have the following decompositions: 



2 DWIA model of (e, e'p) 



da 



- = E + hA , 



(1) 



de' e dn' e dp' 



E = Ka M (v L R L + v T R T + v TL R TL + v TT R TT ) , 



(2) 
(3) 
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Figure 1: Kinematics for the (e, e'p) reaction. The (x, y, z) coordinate system is referred to 
the scattering plane with the z-axis lying along the direction of the momentum transfer q. 
The barycentric system (/, t, n) is referred to the reaction plane: I lies along the direction 
of the ejected nucleon p', the direction n is defined by q x p', and t = n x I. 

om is the Mott cross section, the factor K = mjvp'/(27r/i) 3 , with the nucleon mass, 
and the t> a -coemcients are the usual electron kinematical factors 1 . 

The hadronic dynamics of the process is contained in the exclusive response functions 
R K , which are given as 

R L = R T = W xx + W yy^ ^ 

rtl = ^/2(w Ox + W x0 ), R TT = W yy -W xx , (5) 
RT' = i(w x y - Wv x ), R TU = iV2{W 0y + W y0 ), (6) 

with W^ v the hadronic tensor 

= ^Y,(P^ B \H^\Ar(p's,B\r( q )\A) (7) 

A M B 

constructed from the matrix elements of the electromagnetic nuclear current operator 
J M (q) between the ground state of the target nucleus \A) (assumed to have zero total 

this work we consider the kinematical factors similar to those expressions presented in [48]. Notice 
that these factors differ from the ones of ref. [24] in a global sign for vtl, vtl> and vtt, and an additional 
1 / y/2 factor in the case of the interference TL coefficients 
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angular momentum), and the final hadronic states \p's,B). In what follows we assume 
the residual nucleus to be left in a bound state, hence its wave function can be written 
down in the form \B) = \JbMb) with Jb the total angular momentum. The state \p's) 
represents the asymptotic distorted wave function of the ejected nucleon polarized along 
an arbitrary s-direction, determined by the angles (9 s ,(p s ) referred to the xyz coordinate 
system of Fig. 1. It is given by 

Ip'«> = E ^l 2 \e s As^) | P V), (8) 

u=-l/2 2 

where |pV) is referred to the system with the quantization axis along q and the arguments 
of the rotation matrices are the Euler angles that specify the s-direction. 

Isolating the explicit dependences on the azimuthal angle of the ejected nucleon 0' = 0, 
the hadronic responses can be expressed in the form 



R L = W L (9) 

R T = W T (10) 

R TL = cos(f) W TL + sin <f)W TL (11) 

R TT = cos 20 W TT + sin20 W TT (12) 

R T ' = W T ' (13) 

R TL ' = cos0 l^ TL ' + sin0 W TL ' , (14) 

where the functions W K and W K are totally specified by four kinematical variables, for 
instance {E,u,q,9'}, and the polarization direction {9 S , A0 = — S }. The responses 
with and without tilde refer to their dependence on the spin vector s. As shown below, 



W are purely spin-vector, while W present also a spin-scalar dependence, so only the 
latter survive when the polarization of the ejected nucleon is not measured. 

In the case of (e, e'N) processes, the hadronic response functions are usually given by 
referring the recoil nucleon polarization vector s to the baryocentric system defined by 
the axes (see Fig. 1): / (along the p' direction), n (normal direction to the plane defined 
by q and p', i.e., along q x p') and t (determined by n x I) 2 . It can be shown (see 
Refs. [23, 24] for details) that a total of eighteen response functions enter in the analysis 

— * 

of (e, e'N) reactions. These are given by the decomposition 

W K = \w« npol + W«s n) K = L,T, TL, TT, TL' , (15) 

W K = Wfsi + Wfst, K = TL, TT, T', TL' , (16) 

2 Notice that this notation does not coincide with Refs. [30, 31] where the t direction is denoted as s 
(sideways) 
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where, as mentioned above, only the W^ npol responses survive within the unpolarized 
case. Moreover, W^ pol (referred as fifth response) enters only when the polarization of 
the incident electron is measured. 

Owing to the above decomposition, the response functions (9-14) can be expressed in 
the form R K = R^ npol /2 + H K ■ s, and similarly, the cross section (1-3) can be written as 
a sum of unpolarized and spin-vector dependent terms 

= lz unpol + X.s + h(lA unpol + A.s) (17) 
= hl unpol [1 + P • s + h(A + P' • s)} , (18) 

where the usual polarization asymmetries have been introduced [3]: 



£/ ^-X un po2^ Induced polarization, (19) 
P' = A/ (-T, unpo i) Transferred polarization, (20) 



P 



,2 

A = A unpo i/Y, unpo i Electron analyzing power. (21) 



2.2 Multipole analysis of response functions 

In this section we present the multipole expansion of the response functions to be used in 
our DWIA model. The final expressions, where the sums over third components of angular 
momenta have been performed analytically, result convenient in the present work since 
the computational time can be considerably reduced, specially the calculation concerning 
the MEC. Note that the number of multipoles needed to get convergence increases with 
q, u and up to ~ 36 multipoles are needed for q = 1 GeV/c. The expansion is performed 
following the formalism developed in [49] for exclusive reactions from polarized nuclei. A 
basic difference between the present work and that of ref. [49] lies on the sums performed 
over the third components which are different when initial and/or final state polarizations 
are considered. Here we simply present the final expressions, referring to ref. [49] for details 
on the expansion method and to the Appendix for an outline on the procedure used to 
perform the sum over third components in the present case. 

In order to compute the hadronic tensor in our DWIA model we first perform a 
multipole expansion of the ejected nucleon wave function in partial waves. The final 
hadronic states may then be written 

\p'u,B) = ]T i%* M (f>'){\ulM\ 3 m p )( 3 m p J B M B \J f M f )\{l 3 )J B , J f M f ) , (22) 

IMjrripJfMf 
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where the partial waves (Ij) are coupled to the angular momentum Jb of the residual 
nucleus to give a total angular momentum Jf in the final states |/) = \(Ij)Jb, JfMf). 

The electromagnetic charge and transverse current operators are expanded as sums 
involving Coulomb (C), electric (E) and magnetic (M) tensor operators, 

oo 

p(q) = V^J2 iJ i J Wjo(q) (23) 
j=o 

oo 

J m = -yfa £ i J [J] [Tfm(q) + mfZ 9 (q)} , m = ±1 , (24) 
j=\ 

where, as usual, we assume the transfer momentum q along the z-direction and J m are the 
spherical components of the current operator J. We use the bracket symbol [J] = \/2J + 1 
for angular momenta. Inserting (8,22^23,24) into the hadronic tensor (7), the following 
expansion for the responses W K and W K is obtained 

W K = \w« pol + 2nP 1 1 (cos9 s ) sin(A0)W^, (25) 

for K = L,T, TL, TT, TV, and 

W K = 2-na K [p?(coa9 8 )W£ + Pi(cos^) cos(A0)W£] (26) 

for K = TL, TT, T , TL' , with the coefficient ot K = -1 for TL, TT and a K = 1 for V , TL' , 
and Pj* (cos 6 s ) the Legendre functions. 

The five response functions W^ n ol are the only ones that survive when summing over 
final spins ±s, in which case the 1/2 factor cancels and the unpolarized cross section 
is recovered. The spin dependence is determined from the thirteen reduced response 
functions W^ M (M = 0, 1) and introduced above. Explicit expressions for these 
reduced responses can be written in terms of the reduced matrix elements of the current 
multipole operators 

C CT = ((lj)J B ,J\\Mj\\0) (27) 
E a = ((lj)J B ,J\\ff\\0) (28) 
M a = ((lj)J B ,J\\ifT a9 \\0), (29) 

where we have defined a multiple index a = (l,j, J) corresponding to the quantum num- 
bers of the final states. Note that the initial state \A) = |0) has total angular momentum 
equal to zero, so J/ = J. The response functions involve quadratic products of these 
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c;,c a 


= 


+ Ha'a 


E* a ,E a + M* a ,M a 


= R a'a 


+ ll a'a 


/•:,;-. u, - M*,E a 


- K a'a 




c;,E a 


- K a'a 


4- iJ TL1 




r>TL2 
- K a'a 




E* a ,E a - M* a ,M a 


pTTl 
- K a'a 




E*,M a + M*,E a 


~ U a'a 


, JTT2 



multipole matrix elements which can be decomposed into their real (R^> a ) and imaginary 
(I a' a) Parts: 

(30) 
(31) 
(32) 
(33) 
(34) 
(35) 
(36) 

Expressions for the unpolarized response functions W^ npol in terms of (30-36) are given 
in ref. [41], while the recoil nucleon polarized responses can be written as 

wft = ^E*W*)E p iW**v( n o oW'-<v-< 7 ^ (37) 

W Tl = -J7j2 P J'+Lh{j>Lo(°')l^ P l+l>+Ji®v'v ( ^ \ n ) 
21 J'L a'a \ J 

x + (38) 

~2^E(-i) r+ ^Ui(^)E^W'^f o i -i) 

^ J'L a'a \ / 

( (C/'-z'.j-z-Cv 1 ~~ £.j'-i',j~ilJ'a 2 ) (39) 

~E(-i)^WOE^W^("f f _ L 2 ) 

x (^-iv-i^ 1 - G-iv-i^ 2 ) (40) 

= -^2^E(-i)^UiWE^W'^f o i -i) 

^ J'L a'a \ / 

x^-m.C-G-.^C) (4i) 
wSS = -^2v^E(-ir + ^ ai (^E^W'^( o i 7 _i) 

x (C/'-z'.j-z-Cv 1 ~~ £j'-i',j-ilj'a 2 ) (42) 

wS5 = -^E(-i)^ +L ^L2WE^«W^("f f _ L 2 ) 

x {^.j'-i',.;-^!^ 1 ~ 0'-/',j-/^ct 2 ) (43) 
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L 




J'L o'er \ 

x + 0-v,j-i%?*) ( 44 ) 



W 1M 



^El-il^^ilOE^W^v o i -i 

21 J'L (T'<7 V 



x (^j'-jv-j-^ct'ct 1 ^j'-i',j-i^^f) ? (45) 
where we use the parity functions 

P B ± = (l±(-in/2, ^^(-l) (J '- J)/2 ^, Gj = (-l) (J '" J+1)/2 ^+j (46) 

and the angular dependence of the above responses is determined by the functions h-j JILM (6 
and h-jj, LM (9'), defined through the coupling of two spherical harmonics (see (88)). 

Re[Yj(s)Yjip')] LM = cosM<p' h jj'LM(0')P^(co S 9 s )cos(MA<p) 

M=Q 
J 

+ sinM0' ]T h^ JILM {e')Py {cos 9 s ) sin(MA0) . (47) 

M=0 

Finally, the coefficients $> a > a are derived in the Appendix and are given by eq. (86) selecting 
J = 1. Although the above expressions correspond formally to those denoted as 
and Wjff in [49] for polarized nuclei and J = 1, it is important to point out that the 
coefficients contain the whole information on the polarization distribution of the 

particles. Hence the significance of <& a 'a is clearly different when polarization degrees of 
freedom are considered for the ejected nucleon (present work) or the target nucleus [49]. 

The nuclear structure information in (37-45) is contained in the quadratic forms (30- 
36) of the C, E, M multipoles 3 . Thus the present expansion can be applied to any nuclear 
model of the reaction as far as it provides multipole matrix elements (27,28,29) for high 
enough angular momenta a = (l,j,J). Note that only the responses involving the real 
parts R^, a survive when FSI are neglected since in this case all the C, E, M multipoles are 
strictly real functions. Therefore those responses which depend on the imaginary parts 
are expected to be particularly sensitive to the description of FSI. 

Writing down explicitly the Legendre polynomials involved in the multipole expansion 
(25,26), and comparing with the general expression (15,16), we get the following relation 
between both sets of response functions: 

W? = -2ttW£, K = L,T,TL,TT,TL' (48) 

W t K = 2ti(x k (W£ sin 9' + W£ cos 9') , K = TL, TT, TL', T' (49) 

W t K = 2-na K (W£cos0' - W^sm9') , K = TL, TT, TL', T' (50) 



3 Note that there is a typo in eqs. (40-51) of ref. [49]: the order of J and J' in the three-j should be 
reversed. This error has been corrected in eqs. (37-45) in the present paper. 
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with ax as introduced in (26). 



2.3 Electromagnetic operators and PWIA 

In this work we evaluate the exclusive polarized responses using a semi-relativistic (SR) 
model for describing the electromagnetic one-body (OB) and two-body MEC current 
operators. The OB current has been obtained by a direct Pauli reduction of the fully 
relativistic operator in powers only of the initial nucleon momentum over the nucleon 
mass p/mjv. The dependence on the transfer and final momenta, which can be large [42, 
43, 44, 45], is treated exactly. The SR-OB current in momentum space can be written as 

■Ap',p) = pc + ip so (cos(t) a x -sin0 a y )x (51) 

J X (P,P) = iJmVy + Jc XCOS0 (52) 

J y (p',p) = -iJmOx + Jc X sin 0, (53) 



where x — (p/ m N) sin6>, and (9, 0) are the angles determining the direction of the initial 
momentum p in the (x,y,z) coordinate system. The charge (p c ), spin-orbit (p so ), magne- 
tization (J m ) and convection (J c ) terms shown above include relativistic corrections and 
are given by the following expressions 

Pc = —j=G E , p so = k — -j==== (54) 
V r 2V1 + r 

Jm = VtGm, Jc = Ge , (55) 

K 



where k = q/2rriN, t = \Q 2 \/Am 2 N) and Ge, Gm are the electric and magnetic nucleon 
form factors for which we take the Galster parameterization [50]. 

The two-body MEC operators of pionic (P), seagull or contact (S) and A-isobar kinds, 
displayed in the Feynman diagrams of Fig. 2, have been also obtained by making use of a 
SR approach leading to simple prescriptions that include relativistic corrections through 
a multiplicative factor (see Refs. [40, 45, 46] for details on the SR expansion method) 



jMEC _ _} -iMEC /r R \ 

J SR — / -, ; NR 5 \ 00 ) 
Vl + T 



where Jjvif is the traditional non-relativistic MEC operator. 

The expressions for the reduced matrix elements of the OB and MEC multipole oper- 
ators (27,28,29) in the shell model are given in Refs. [10, 11, 43] except for the relativistic 
correction factors appearing within the SR operators. The somewhat complex structure 
displayed by these multipoles makes it not possible to predict the relative importance of 
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(a) (b) (c) (d) 



Figure 2: MEC diagrams contributing to the two-body current of this work: seagull (a), 
pion in flight (b) and A (c,d) currents are considered including relativistic corrections. 

each contribution separately without explicit numerical evaluation, even in the case of the 
OB current. 

Although in this work we perform a DWIA analysis of the response functions, we 
may take advantage of the significant simplifications introduced within the plane wave 
impulse approximation (PWIA), where analytical expressions for the response functions 
can be obtained [30, 31]. First, for intermediate to high values of q, the PWIA approach 
is expected to provide reasonable results, thus the analytical PWIA expressions allow 
us to estimate the contributions of the different pieces of the currents to the polarized 
response functions. Second, since the PWIA results should be recovered using the present 
multipole expansion in the limit of no FSI, the comparison between our calculation and 
the analytical PWIA responses makes it possible to fix the number of multipoles needed 
to get convergence. 

Hence, within PWIA, the matrix element of the OB current is written as 

(p's,B\J»(ci)\A) = ^V;, 1 _(s)J»(p , ,p)p, p (B\a p ,p\A) , (57) 

where a Pt/ 3 is the annihilation operator corresponding to a particle with momentum p 
and spin projection /3 referred to the quantization axis. Inserting this expression into the 
hadronic tensor (7), and following the procedure described in [43], we obtain 

W" = ^m N p'w^(p', p, s)M s (p) , (58) 

where we have defined the polarized single-nucleon tensor 

w^(p',p,s) = £ P^( S )r(p',p)^^(p',p)*, a P a ,i( S ). (59) 

aa'P' 
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In the case of interest here, a closed-shell nucleus, the scalar momentum distribution 
M 5 (p) for nucleon knock-out from a shell nlj is given by 



M*(p) = 2J ^R 2 (P) (60) 

with R(p) the radial wave function of the hole in momentum space. 

Using the current matrix elements (51,52,53) one can compute in the factorized ap- 
proximation (58) the response functions (4-6). From these results the PWIA expressions 
for the reduced response functions can be identified. Expressions for the unpolarized re- 
sponses in PWIA were given in [43, 49]. In the case of the polarized responses, from the 
total of eighteen, only five survive in PWIA. These are given by 



< = ^2J c J mX M s (p) (61) 

< = ^2£M*(p) (62) 
W^' = ^2V2(p c J m - Pso J cX 2 )M s (p) (63) 
W™ = ^-2y/2p e J m M 8 (p) (64) 

= ^-2V2 Pso J mX M s ( P ), (65) 



where we have used the factors introduced in (54,55). Notice that all the L,T,TL and 
TT-type polarized responses are zero in this approximation. 

3 Results 

In this section we present results for selected recoil nucleon polarization observables cor- 
responding to proton knockout from the pi/2 and p^/2 shells in 16 0. In particular, we 
restrict ourselves to the analysis of all the polarized response functions, including the fifth 
response Wq L that does not depend on the nucleon polarization and only enters when 
the initial electron beam is polarized, and the transferred polarization asymmetries P( t n - 
The study of cross sections and induced polarizations will be presented in a forthcom- 
ing publication [51]. Two different kinematical situations corresponding to (q, u;)-constant 
kinematics (also referred as quasiperpendicular kinematics) have been selected: i) q — 460 
MeV/c, uj = 100 MeV, and ii) q = 1 GeV/c, u = 450 MeV. In both cases the value of the 
transfer energy uj corresponds almost to the quasielastic peak. 

In this work our main interest is focused on the role of the two-body MEC operators 
upon the recoil nucleon polarization observables, trying to identify kinematical conditions 
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for which these effects can be important; however, a brief excursion on the FSI effects 
is also presented. All the calculations have been done within the formalism described in 
the previous section, i.e., semirelativistic expressions for the one- and two-body current 
operators and a multipole expansion method have been used. The number of multipoles 
needed has been fixed by comparing the DWIA results, in the particular case of no FSI, 
with the exact factorized PWIA responses (61-65). Convergence in the multipole analysis 
is obtained with J max = 30 for q = 460 MeV/c, and J max = 35 for q = 1000 MeV/c. 
Finally, in all of the results which follow, the kinematics of the ejected nucleon is treated 
exactly. 

3.1 Polarized response functions 

Here we analyze the thirteen responses defined in (48-50) which arise from the ejected 
nucleon polarization, plus the "fifth" response function Wj L . Results are displayed in 
Figs. 3-12. A similar analysis for the unpolarized responses L, T, TL and TT has been 
performed recently in [40]. 

3.1.1 Effects of FSI 

We start our discussion with the effects of FSI on the polarized responses. A study of 
the dependence of the unpolarized responses on the particular FSI model was already 
presented in [38]. In Fig. 3 we show the eight induced polarized responses for proton 
knock-out from the lpi/2 shell in 16 as function of the missing momentum p. Kinematics 
corresponds to q = 460 MeV/c and oo = 100 MeV. The five transferred polarized responses 
plus the fifth one (T" and TV types) are displayed in Fig. 4. Similar results are obtained 
for the IP3/2 shell and thus they are not shown here. In all of these results we use bound 
wave functions obtained as solutions of the Shrodinger equation with a Woods-Saxon 
potential, with parameters taken from ref. [42]. For the final states we use solutions for 
two different optical potentials. Solid lines correspond to calculations performed with the 
Comfort and Karp potential [52], which was originally fitted to elastic proton scattering 
from 12 C for energies below 183 MeV. We have extended it to 16 O by introducing a 
dependence A 1 / 3 in the radius parameters. The results shown with dashed lines have 
been computed with the Schwandt potential [53], which also has been extrapolated here 
for 16 O since it was originally fitted to higher mass nuclei. 

The induced polarized L, T, TL, TT and the fifth response functions, which are zero in 
absence of FSI, are expected to be highly sensitive to the details of the particular optical 
potential considered, and in particular, to the spin-orbit term in the potential. In this 
sense notice the significant difference introduced by both potentials in the case of the 
polarized TL and TT responses (fig. 3), while the FSI discrepancy gets smaller for the 
fifth response function (fig. 4) and is considerably reduced for and ■ 
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p [MeV/c] p [MeV/c] 

Figure 3: Induced polarized response functions (L, T, TL and TT-type responses) for 
proton knock-out from the lpi/2 shell in 16 O, with momentum transfer q = 460 MeV/c 
and energy transfer uo = 100 MeV. Solid lines are the DWIA results using the optical 
potential of Comfort and Karp; dashed lines have been computed with the Schwandt 
optical potential. 
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Figure 4: The same as fig. 3 for the fifth and transferred polarized response functions (T" 
and TL'-type responses). With dotted lines we show also for reference the PWIA results. 
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The five transferred polarized responses which survive in PWIA (61-65), depend less 
on the details of the potential, being mostly affected by the central imaginary part of it. As 
known, these responses enter in the case in which also the initial electron is polarized, and 
they contribute to the transferred nucleon polarization asymmetry We observe (fig. 4) 
that both potentials lead to close results, differing by less than 10% for the dominant 
responses W^ L> , Wf' and Wf L ', while the largest differences are shown for Wf\ which 
is however very small. Similar results are found for the P3/2 shell. The sensitivity shown 
by some polarized responses to the details of the potential, makes these observables of 
special interest to disentangle between the different models of FSI that can fit reasonably 
well the unpolarized cross sections. 

3.1.2 Effects of MEC 

The impact of MEC on the recoil nucleon polarized responses is shown in Figs. 5-12. 
In each panel we compare the distorted wave responses evaluated by using the OB cur- 
rent only (dotted line) with the results obtained when including also the two-body MEC 
operators considered in fig. 2, namely, the seagull or contact (OB+S) current (dashed 
lines), the contact and pion-in-flight (OB+S+P) currents (dot-dashed lines), and finally, 
including also the A current (solid lines), denoted as (OB+MEC). Results in figs. 5-8 cor- 
respond to kinematics q = 460 MeV/c, uo = 100 MeV (kinematics I), whereas in figs. 9-12 
we present the responses evaluated at q — 1 GeV/c, uo = 450 MeV (kinematics II). For 
both kinematics proton knock-out from the pi/ 2 (figs. 5,6 and 9,10) and ps/ 2 (figs. 7,8 and 
11,12) have been considered. 

Let us discuss first the results for kinematics I (figs. 5-8). Here we observe that the 
global sign of the polarized T, TL and TT responses changes when comparing the p±/ 2 
(fig. 5) andp 3 / 2 (fig. 7) shells. The same occurs for kinematics II. Concerning MEC effects, 
the various polarized responses display different sensitivities to the two-body component 
of the nuclear current. Apart from the pure longitudinal response W£, which shows no 
dependence on MEC because the "semi-relativistic" MEC expressions only include the 
leading transverse components, the role of MEC on W£ is shown to be similar to the 
one found for the unpolarized T- response in [40]: the enhancement (in absolute value) 
produced by the S current is partially cancelled by the reduction introduced by the P 
current; the A current gives rise to an additional reduction, leading to a global decrease 
of the response of the order of ~ 10% at the maximum. This effect being similar for 
both shells (figs. 5 and 7). 

Larger MEC effects are found for some of the induced polarized TL responses, par- 
ticularly for Wj L where the A current produces a very significant modification of the 
response, changing even its shape in the region close to p ~ 100 MeV/c. Note that, al- 
though the global effect introduced by the A in this response is similar for both shells, in 
the case of the pi/ 2 there is a large increase, whereas for p 3 / 2 the response is significantly 
reduced in absolute value. It is also interesting to point out that the A current plays also 
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O, lpi/2, q = 460 MeV/c, u = 100 MeV 




Figure 5: Induced polarized response functions for proton knock-out from the lpi/2 shell in 
16 0, with momentum transfer q = 460 MeV/c and energy transfer u = 100 MeV. Dotted 
lines are the DWIA results using only the OB current operator; dashed lines include in 
addition the seagull current (OB+S); dot-dashed include the seagull plus pionic currents 
(OB+S+P); finally solid lines represent the total result (OB+MEC) including also the A 
current. 
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Figure 6: Fifth and transferred polarized response functions for 16 0. Kinematics corre- 
sponds to momentum transfer q = 460 MeV/c and energy transfer u = 100 MeV. The 
meaning of the lines is the same as in Figure 5. 
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the most important role for the response, this being clearly shown in the case of the 
p 3/2 shell. 

The role of MEC on the three polarized TT responses shows a very different behaviour 
for the two shells considered. In the case of the pi/2 (fig. 5), the global effect of MEC 
is a very significant reduction of the responses, particularly for W^ T (~ 20%) and W^ T 
(~ 30%), being the separate contributions of the S, P and A currents of rather similar 
importance. Note that the contributions introduced by the S and P currents have opposite 
signs for the pi/ 2 and p^/2 shells. As a consequence, for the p^/2 shell (fig. 7) the large 
enhancement (in absolute value) produced by the S current is almost cancelled exactly 
by the contributions of the P and A currents, so the net MEC effect is almost negligible 
for the three TT responses. 

The transferred polarized responses (T" and TL'-type responses) are shown in Figs. 6 
and 8. From these results we find in general a small effect of MEC, less than ~ 5%. 
An exception is Wj' where the role of A gives rise to an important reduction of the 
response; however notice that Wj is very small, of the order of ~ 10% compared with 
Wf and Wj L , and hence difficult to measure. The anomalous smallness of the response 
Wj was already discussed in detail in [31] within the context of the PWIA and differ- 
ent non-relativistic reduction schemes. This result can be also understood within the 
multipole analysis performed in this work by taking into account the general relations 
given in eqs. (49,50) and the explicit expressions obtained for the multipole functions in 
PWIA (61,62). Since we are close to the quasielastic-peak, the angle 9' is close to zero for 
moderate missing momentum, so the biggest contribution comes from the factor multi- 
plied by cos 9' in eqs. (49,50). This factor is = 0(1) in the case of the /-responses, and 
Wf[ = 0(x) f° r the t-response. Precise values of the nucleon form factors and kinematical 
variables can be introduced in these equations to verify the exact relation between the / 
and t components in PWIA, which is not very different from the distorted-wave results 
of figs. 6 and 8. 

Results for higher momentum and energy transfer, q = 1000 MeV/c and uo = 450 
MeV (kinematics II), are shown in Figs. 9-12 for the two p shells in 16 0. This kinematics 
corresponds to the experimental setting of [16, 22] where Q 2 = —0.8 (GeV/c) 2 . Obviously 
in this case relativity is expected to play a more important role and in fact, studies 
within the relativistic distorted wave impulse approximation (RDWIA) [13] have proved 
the importance of these effects. The present SR model, although lacking some of the 
relativistic ingredients inherent in the RDWIA, incorporates exact relativistic kinematics 
for the ejected nucleon, a SR expansion of the current which can be used for high q 
values, and finally, the use of the Schrodinger-equivalent form of the S-V Dirac global 
optical potential of [54], including the Darwing term in the wave function. The validity 
of the expansion procedure used in the SR model was tested in [13] where unpolarized 
observables evaluated within the SR approach were compared with a RDWIA calculation 
for this kinematics. 

The discussion of the results presented in figs. 9-12 follows similar trends to the ones 
already presented for kinematics I, so here we simply summarize those aspects which can 
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Figure 7: The same as Fig. 5 for proton knock-out from the 1^3/2 shell. 
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16 0, lp 3/2 , q = 460 MeV/c, u = 100 MeV 
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Figure 8: The same as Fig. 6 for proton knock-out from the 1^3/2 shell. 
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O, lpi/2, q = 1000 MeV/c, oo = 450 MeV 




i i i 

100 200 300 400 
p [MeV/c] 



Figure 9: The same as Fig. 5 for q 
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1000 MeV/c and uj = 450 MeV. 
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O, lpi/2, q = 1000 MeV/c, u = 450 MeV 
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Figure 10: The same as Fig. 6 for q = 1000 MeV/c and u = 450 MeV. 



23 



16 0, lp 3 /2, q = 1000 MeV/c, u = 450 MeV 




p [MeV/c] p [MeV/c] 

Figure 11: The same as Fig. 5 for the p 3/2 shell, q = 1000 GeV/c, and uo = 450 MeV. 
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16 0, lp 3/2 , q = 1000 MeV/c, u = 450 MeV 




p [MeV/c] p [MeV/c] 

Figure 12: The same as Fig. 6 for the p 3/2 shell, q = 1000 MeV/c and u = 450 MeV. 
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be of more relevance. As shown in figs. 9-12, the general effect introduced by MEC is 
a global reduction of the responses (in absolute value) whose magnitude depends on the 
specific response, being of the order of a few percent for W^ t L n and W^ L tn , larger for 
and Wj" t (particularly because of the A-contribution) and the largest for Wf T : 
where the reduction (basically due to A) is about ~ 20 — 25%. Note however that the 
response Wf T is the smallest one and so hardly measurable. 

The dependence of MEC effects on the momentum transfer shown in the results of 
figs. 5-12 is consistent with the findings of ref. [40] for the unpolarized responses. In 
general the importance of MEC decreases with q. This is in accord with the results for the 
T response in the lp— lh channel in the case of quasielastic inclusive (e, e') reactions [45, 
46]. This behaviour can be roughly understood from the relativistic expressions for the 
particle-hole transverse current matrix elements J T (p',p) in Fermi gas [45, 46], and also 
from the traditional non relativistic expressions. At the non relativistic level, the OB 
current is dominated by the magnetization contribution which goes as ~ q. On the 
contrary, MEC present a much more complex dependence on q and on the momenta 
of the two holes involved: the missing momentum p and an intermediate momentum 
k which should be integrated. Moreover, MEC also contain pion propagators involving 
inverse squared pion momenta. For high q, a crude estimation of the (transverse) seagull 
and pion in flight currents is shown to behave as ~ q/(q 2 + m 2 ), w hile the A current 
goes as ~ q 3 /{q 2 + m 2 ), hence the latter clearly dominates, which is in accord with the 
results shown here. Once the nN form factor, which becomes smaller when high momenta 
are probed, is added to the two-body currents, we find the OB contribution to dominate 
over the MEC. At the relativistic level the above dependences on q change. In [44] it 
was demonstrated that if the form factors are neglected, then the OB, seagull and pionic 
currents grow asymptotically as y/q. Thus the inclusion of tcN form factors is essential 
for the dominance of the OB current This conclusion however applies to the response 
functions only for low missing momentum, since for other observables such as the Atl 
asymmetry [40] and the polarization asymmetries (see below) larger effects are found for 
high values of q and missing momentum. 

3.2 Transferred polarization asymmetries 

Apart from the response functions, other observables of special interest are the nucleon 
polarization asymmetries introduced in eqs. (18-21). These observables are given as ratios 
between polarized and unpolarized responses, where one hopes to gain different insight 
into the underlying physics from what is revealed through the responses themselves. As 
already mentioned in the introduction and in order to clarify the discussion, here we 
restrict ourselves to the analysis of the transferred polarization asymmetries P( t n , which 
only enter with polarized incident electrons and persist in PWIA. Induced polarization 
ratios P^t, n — that do not depend on the polarization of the incident electron and are 
zero within the plane wave approach — , and total cross sections will be analyzed in a 
forthcoming publication [51]. 
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Following the discussion presented for the responses, here we first study the effects 
introduced by FSI and later on we focus on the role of MEC. 



3.2.1 Effects of FSI 



In Figs. 13-14 we present the results obtained for the transferred polarization asymmetries 
corresponding to proton knockout from the p\/2 and P3/2 shells in 16 O. Kinematics has 
been selected as (1), i.e., q = 460 MeV/c and uj = 100 MeV. Results for kinematics (II) 
follow the same general trends although FSI effects are in general less important because 
of the higher momentum transfer involved. The longitudinal P{ and transverse (sideways) 
P[ components are shown in fig. 13 for electron scattering angle fixed to 9 e = 30° (forward 
scattering) and three values of the proton azimuthal angle = 0, 90° and 180°, while the 
normal polarization P' n is displayed in fig. 14 for <fi = 90° (notice that P' n is zero for 
co-planar kinematics). Although not shown here for brevity, we have also explored the 
behaviour of the transferred polarization ratios at backward scattering angle (9 e = 150°). 
As known, the purely transverse responses dominate at backward angles, whereas all of 
the kinematical factors that enter in the description of (e, e'p) reaction are of similar order 
at forward angles. In [30] forward scattering angles were proved to enhance significantly 
the sensitivity to dynamical relativistic effects. Concerning FSI and MEC, the discussion 
of the results for 9 e = 150° follow similar trends to the ones presented here for 9 e = 30°. 

The PWIA calculation (dotted line) is compared with DWIA results using the two op- 
tical potentials already presented in the previous section, i.e., Comfort & Karp (solid lines) 
and Schwandt (dashed lines). First, notice the difference between PWIA and DWIA re- 
sults. Within the plane wave approach, the responses factorize and hence the polarization 
ratios depend only on the single-nucleon responses, being cancelled the whole dependence 
with the momentum distribution. This means that PWIA results are identical for the two 
p-shells considered. Moreover, polarization ratios in PWIA may be written in the general 
form 

D , Oi + hx + 0(x 2 ) , . . , (aa , 

p - = e, + d , x + o ix iy fOTI = '•«•"• (66) 



where x — p/ m N sin 9, already introduced in (51-53), is the parameter in the SR expansion 
of the nuclear current. For low missing momentum, the above fraction has a linear 
dependence on x pl us a small correction of order \ 2 which breaks linearity for higher 
p. 

For low missing momentum values p < 200 MeV/c, the effects introduced by FSI are 
small, being almost negligible at the maximum of the momentum distribution (p ps 100 
MeV/c). This result is expected because of the global reduction of the polarized response 
functions produced by FSI: of the order of ~ 30% (fig. 4). This is somewhat similar to the 
behaviour shown by the unpolarized responses [38]. Hence, although not exact because 
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p [MeV/c] 

Figure 13: Transferred polarization ratios P[ and P[ for proton knock-out from the p shells 
in 16 0, and q = 460 MeV/c, u = 100 MeV. The electron scattering angle is 6 e = 30°, and 
results are shown for three values of the proton azimuthal angle = 0, 90°, 180°. Solid 
lines have been computed in DWIA with the Comfort and Karp potential, dashed lines 
with the Schwandt potential, and finally the dotted lines are the PWIA results. 
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of the slightly different sensitivities to FSI shown by the various responses, a kind of 
cancellation of FSI between the numerator and denominator in the polarization ratios 
occurs for low p. From results in fig. 13, one also observes that FSI effects are slightly 
bigger in the case of the P3/2 shell, particularly for P{ and <p = 180°. The reason of this is 
connected to the much less reduction that FSI cause upon the unpolarized TL response 
forp 3 / 2 (see [38] for details). 

For high missing momentum the DWIA polarizations deviate significantly from the 
PWIA results, showing a very pronounced oscillatory behavior which may even give rise 
to a change of sign in the polarizations. This is a clear indication that for high momentum 
the effects of FSI are not simply a global reduction of the responses due to the imaginary 
part of the potential, but on the contrary, each response turns out to present a peculiar 
sensitivity to the interaction. As shown in figs. 3 and 4, this is hardly visible in the separate 
response functions because of the smallness of the momentum distribution for high p. It is 
important to point out that the oscillatory behaviour presented by the polarization ratios 
is a direct consequence of the breaking of factorization property. This issue was already 
studied at the level of the plane wave approach taking care of the dynamical relativistic 
effects introduced by the lower components of the bound Dirac spinors [30]. A general 
analysis of factorization within the context of the RDWIA and different non-relativistic 
approximations is presently in progress [32]. 

Focusing on the results presented in fig. 13, we observe that the shape and magnitude 
of both polarization asymmetries, P[ and P/, are similar for the two p-shells. In the 
particular case of <fi = 90° (out-of-plane kinematics) the ratio P' t is very small, almost 
negligible for low missing momentum. This is expected since only the response Wj , 
which is very small, contributes to P[ in that situation. For co-planar kinematics a large 
discrepancy between the results obtained at <fi = 0° and <fi = 180° exists. As shown 
by eqs. (13-16), the numerator in the ratios P/, % = l,t, is given through the linear 
combination vx'Wf + vtl'WJ l cos0, with the kinematical factors being vt> = 0.27 
and v T l' = —0.18 for the kinematics considered here (I). Hence, from the transferred 
polarization asymmetries measured at = 0° and = 180°, the separate responses Wf' 
and Wj L ' could be extracted. 

Comparing the solid and dashed lines in fig. 13 we conclude that the uncertainties 
introduced by the optical potentials selected are rather small. For low momentum transfer 
these differences are negligible in contrast to fig. 4 where some responses are shown to be 
affected appreciably by the optical potential. This again is an outcome of the fact that the 
differences between the responses computed with these potentials are of the same size in 
numerator and denominator and they tend to cancel when taking the quotient to compute 
the polarizations. Both sets of results start to differ for p > 300 MeV/c. Note however 
that for high p-values other relativistic effects coming from the dynamical enhancement 
of the lower components in the wave functions, not included in the present model, may 
also contribute significantly to the oscillatory behaviour of the polarizations [30, 32]. 

The case of the normal polarization transfer P' n (fig. 14), present some peculiarities not 
observed for P( t . First the difference between PWIA and DWIA results is rather constant 
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Figure 14: The same as Fig. 13 for the polarization transfer component P' n and <fi = 90°. 

for the two shells in the whole range of missing momentum. In particular, the distorted 
wave approach leads to results which are very similar to the ones obtained within PWIA 
in the case of the ps/2- In addition, the strong oscillatory behaviour due to FSI and shown 
for P( t (fig. 13) does not appear here, being the differences introduced by both optical 
potentials small. These results could promote this observable P' n , that can be obtained in 
out-of-plane experiments, as a good candidate in order to study properties of the reaction 
without being much affected by FSI. 

3.2.2 Effects of MEC 

In Figs. 15-17 we present the effects introduced by MEC upon the transferred proton 
polarization components for the two p shells in 16 O, and for the two kinematics considered 
above. In the case of kinematics (I), i.e., q = 460 MeV/c (fig. 15), where we use the optical 
potential parameterized by Comfort & Karp, MEC effects are shown to be small in the 
whole missing momentum range and similar for both shells. Note also that the role played 
by MEC is of the same order of magnitude or even smaller than the uncertainty introduced 
by the use of the two optical potentials (fig. 13). The smallness of MEC effects on the 
polarization asymmetries comes from their effective cancellation when taking ratios of 
response functions. 
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Results for higher momentum transfer, q = 1000 MeV/c (kinematics II), are shown 
in Fig. 16 for the Schrodinger-equivalent form of the S-V Dirac global optical potential 
of ref. [54]. As in the previous case, MEC effects are small for low missing momentum; 
however they tend to increase significantly for higher p-values due to the A exchange 
current, inducing a softening of the transferred polarization asymmetry, which makes its 
oscillatory behaviour to appear at slightly lower momenta (see for instance the important 
MEC effects observed for P[ at <fi = 0, particularly in the region close to the minimum 
p ~ 300 MeV/c). The present results indicate that the response functions entering into 
the polarization ratios are importantly affected by MEC, mainly due to the A current, for 
high momentum transfer. Large effects of this kind have also been found recently for the 
Atl asymmetry obtained from the analysis of unpolarized (e,e'p) reactions corresponding 
to the same kinematics II [40]. 

The normal polarization P' n is shown in Fig. 17 for the two kinematics and both 
shells. MEC effects follow the same trends as those observed for P{ and P' t : they increase 
significantly for high momentum transfer (q = 1000 MeV/c) and high missing momentum. 
In contrast with the P[ and P[ cases, here the relative contributions of the separate MEC 
currents depend on the specific kinematics and p-shell selected, playing the seagull and 
pion-in-flight currents an important role. 

To finish we present in Fig. 18 the asymmetries P{ and P[ evaluated for q = 1000 
MeV/c, uj = 450 MeV and electron incident energy e e = 2450 MeV. This kinematics 
corresponds to a recent experiment performed at TJlab. We compare our calculations 
with the experimental data presented in [22]. The azimuthal angle in this experiment 
was = 180°. Note the change of sign of P[ with respect to the results of Fig. 16, 
due to the opposite definitions of the normal plane (and hence of the t component) for 
4> = 180° (in our case the normal plane for = 180° would point down in Fig. 1, while 
in ref. [22] it was chosen along the up direction). Results for the p and s shells in 16 O 
are shown from left to right. Although being aware of the possible modifications that the 
"dynamical" relativistic ingredients [13, 27, 22] may introduce in the present calculations, 
we are rather confident that the results in fig. 18 give us a clear indication of how much 
the DWIA calculation is expected to be modified after including the two-body (MEC) 
contributions (compare dotted with solid lines). As noted, whereas the contribution of 
MEC over P' t is negligible, they give rise to a slight reduction of P(, which is well inside 
the experimental error except for the S1/2 shell. Comparing the results for the two p-shells 
we observe that our model describes better the case of the ps/2- This is in agreement with 
the findings in [13, 40] concerning the A TL asymmetry. The particular case of s 1 / 2 shows 
that the experimental data for P[ are well reproduced by the calculations, which however 
understimate very significantly the data for P[. 

In order to clarify the importance of FSI, in fig. 18 we also show with dot-dashed lines 
the results corresponding to the global OB+MEC calculation but without including the 
spin-orbit term of the optical potential, i.e., using a phenomenological optical potential 
consisting only of a central part. As shown, the corresponding polarizations present some 
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Figure 15: MEC effects over the transferred polarization asymmetries P{ and P' t for 
proton knock-out from the p shells in 16 0, and q = 460 MeV/c, u = 100 MeV. The 
electron scattering angle is 9 e = 30°, and results are shown for three values of the proton 
azimuthal angle = 0, 90°, 180°. The meaning of the lines is the same as in Fig. 5. 
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Figure 16: The same as Fig. 15 for q = 1000 MeV/c and u = 450 MeV. 
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Figure 17: The same as Figs. 15 and 16 for the normal polarization transfer component 
P' n and <p = 90°. 
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kind of "linearity", being similar to the PWIA results. This is expected since the spin- 
orbit interaction is the main responsible of the oscillatory behaviour of the polarization 
ratios, apart from the "dynamical" relativistic effects. 

The ratio P//P/, shown in the bottom panels of fig. 18, has been proposed as a suitable 
observable for getting information on nucleon properties inside the nuclear medium [55]. 
From inspection of fig. 18 we find that MEC produce a tiny reduction of this observable, 
particularly for low missing momentum, being larger as p goes up. 

All of the above results have been computed using the Galster parameterization for 
the nucleon form factor. It is of interest to know the dependence of our results with the 
nucleon structure, hence we have also calculated the OB+MEC polarization asymmetries 
assuming the Gari-Krumplemann (GK) form factor parameterization [56]. The results are 
shown with dashed lines. The GK parameterization was used by Udias et al. [57] within 
the context of the relativistic calculations presented in ref. [22]. The P[ computed with 
GK form factors is increased with respect to the solid lines, being a little bit closer to the 
experimental data. Let us remind that the effects of MEC lead to a global reduction of 
all of these polarization observables, hence the OB calculation using the GK form factors 
would be clearly located above the corresponding results including MEC (dashed lines). 
This makes our present results to come closer to the relativistic ones of [22]. Note also 
that the uncertainty introduced by the nucleon form factor parameterization shows up in 
P(, being negligible for P' t . 

To finish the discussion, it is also interesting to point out that the behaviour shown 
by the P' t data, growing with p for pi/ 2 and the reverse for p 3 / 2 , does not agree with the 
theoretical results which increase with p for both shells. This is in accordance with other 
relativistic calculations [22]. For p = 140 MeV/c our predictions for P[ in the case of the 
P1/2 shell clearly underestimate the data; as already mentioned, other relativistic effects 
coming from the lower components of the Dirac wave functions, not considered here, may 
also play a significant role. 

3.3 Comparison with previous works 

Concerning previous calculations of MEC in (e, e'p) reactions, in [38, 40] comparisons for 
unpolarized observables obtained with the present model with those of Refs. [34, 35, 37, 39] 
were presented. Next we summarize the differences of MEC effects on recoil polarization 
observables between the present work and Refs. [33, 36]. 

1. Bofli and collaborators [33] find for intermediate q values large MEC effects on the 
separate polarized responses (reduction of the order of 15-30% or even larger), being 
the A current the main contribution. We get in general smaller and qualitatively 
different effects for this kinematics, the seagull and pion in flight currents being in 
our case as important as the A. Concerning the transfer polarization ratios, they 
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Figure 18: Transferred polarization asymmetries P( and P/, and their quotient P[jP{ for 
proton knock-out from the three shells in 16 for q = 1000 MeV/c and uo = 450 MeV. 
The incident electron energy is e e = 2450 MeV and the proton azimuthal angle = 180°. 
The meaning of the lines is the following: dotted are the DWIA calculation with OB 
current only; solid are the total OB+MEC result; dot-dashed are also the OB+MEC 
result but without the spin-orbit term of the optical potential. These three curves have 
been obtained using the electromagnetic nucleon form factors of Galster. Finally, dashed 
lines are the total OB+MEC result using instead the form factors of Gari-Krumplemann. 
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find P( to be the most sensitive one, with a 20% decrease due to the A for low 
missing momentum p < 200 MeV/c. In our calculation MEC effects are clearly less 
important for these missing momentum values. 

2. Ryckebusch et al. [36] do not present the separate l,t,n response functions. In 
general they find small MEC effects, as we do, in the transferred polarizations for 
low p < 300 MeV/c. These effects being larger as q and p increase. Comparing 
specifically our results to theirs for kinematics II, we observe that the OB results 
clearly differ due to the different treatment of FSI, while somewhat larger and 
qualitatively different MEC effects are found in this work. 

Be it as it may, since the different treatment of FSI and of the current operators in [33, 
36] and in the present work produces discrepancies already at the impulse approximation, 
it is hard to draw general conclusions on MEC effects beyond the fact that in [33] MEC 
lead to excesively large contributions compared with us, while their small size in [36] is 
in accord with our calculation. 

4 Summary and conclusions 

In this paper we have presented a distorted wave model of (e, e'p) reactions which goes be- 
yond the impulse approximation with the inclusion of two-body meson exchange currents. 
Relativistic kinematics to relate the energy and momentum of the ejected proton is used 
and the currents are derived through an expansion in powers of the missing momentum 
over the nucleon mass. Explicit expressions of the polarized response functions in a general 
multipole expansion method are given. Results for the responses and transferred polar- 
ization asymmetries have been obtained for proton knock-out from the different shells in 
16 for quasiperpendicular kinematics with the transfer momentum fixed to q — 460 and 
1000 MeV/c. FSI have been considered in each case by using different optical potentials. 

One of our primary goals has been to estimate, within our present model, the validity of 
the impulse approximation by analyzing the effect of MEC on the different recoil nucleon 
polarized observables. Thus we compare the standard DWIA results, obtained using 
only the OB current, with the "full" calculation which includes the MEC. We have also 
explored the role played by the particular description of the FSI, hence we compare the 
results obtained by using different optical potentials which have been widely considered 
in the literature: Schwandt [53] and Comfort & Karp [52] parameterizations. For higher 
energy we have used instead the Schrodinger-equivalent form of a Dirac optical potential. 

From our present studies we may summarize and conclude the following: 

1. The induced T, TL and TT polarized responses are particularly sensitive to the 
details of the optical potential, allowing them, specially the TT ones, to constrain 
the theoretical model for FSI. The transferred polarized responses (T" and TL'), 
which survive in PWIA, show a much less sensitivity to the interaction. 
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2. In general, MEC effects over the transferred T', TV polarized responses for q = 460 
MeV/c and moderate missing momentum (p < 300 MeV/c) are rather small (< 5%), 
and tend to increase as q goes higher, being of the order of a ~ 10% reduction (due 
mainly to the A current) in the particular case of Wj' and q = 1 GeV/c. 

The role of MEC gets clearly more important for the induced T, TL and TT polar- 
ized responses. Emphasis should be placed on Wf T and W^ T which are reduced at 
the maximum by ~ 20% and ~ 30%, respectively, for q = 460 MeV/c and for the 
P1/2 shell; notice however that these effects are negligible in the case of ps/2- For 
q — 1 GeV/c the role of MEC diminishes. 

3. FSI give rise to an important deviation of the transferred polarization asymmetries 
P'l and P[ with respect to the PWIA results, showing a very pronounced oscillatory 
behaviour that starts for p > 200 MeV/c. This behaviour does not appear in the 
component P' n . The uncertainties introduced by the optical potentials are rather 
small for the missing momentum region analyzed. 

4. MEC effects on P[ and P( are negligible for q = 460 MeV/c and increase for q — 1 
GeV/c, especially for p > 200 MeV/c. The role of MEC on P' n is clearly more 
important. 

Finally we are confident that the significant sensitivity shown by some polarized ob- 
servables to MEC, particularly to the A current, will be maintained within the scheme 
of a "fully" relativistic calculation which takes care of relativistic ingredients, such as the 
dynamical enhancement of lower components, not included in the present model. Work 
along this line is in progress. 
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A Sum over third components and reduced response 
functions 

In this appendix we perform the sum over third components in the multipole expansion 
of response functions and give their explicit expressions in terms of the reduced matrix 
elements of the multipole operators. 
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First we write the response functions (4-6) in terms of the spherical components of 
the current matrix elements J ±1 = (p's, B\ J±\\A) using the hadronic tensor (7) 

R L = ^EP*P (67) 

R T = ^EO^-ir + l^il 2 ) (68) 
R TL = - 2 l Re ]>>*(J +1 - J_ x ) (69) 



R 



TT 



K 

K 

■)TL' - 1 



= ^E^V+i + ^V-i) (70) 



i? TL = -2-ReEp*U+i + ^i) (71) 



it> 2 



^E(l^i| 2 -I^i| 2 ) • (72) 

(73) 



is: 



Inserting the multipole expansion for the charge and current components as given in (23,24) 
we find that each response can be written as a sum of terms of the type 

B^r^^Y,(P s ^BM B \f\j lm i q )\A}^p's,J B M B \f Jm (ci)\A}, (74) 

A M B 

where Tj, m , and Tj m represent in general the Coulomb, electric or magnetic multipole 
operators. Introducing now the multipole expansion (22) corresponding to the final state, 
which is polarized along an arbitrary direction s (8), we get 

b?!7 = ^EE^f(-)^ 2) (-)* 



M B v'v 

X 



E /^v(pO(^W|/m;)(/m;j B M B |jW)((zy)j B ,jW|T> m ,|A)* 

l'M'j'm' p A 

x E r^ M (pO(^/M|jm p )(jm p J B M B |Jm)((/jVi?,^l^rn|^), (75) 

IMjrrip 

where we have used Jf — J and Mf = m, since the initial nucleus has total angular 
momentum equal zero. 

Using the Wigner-Eckart theorem for the matrix elements of tensor operators between 
states of definite angular momenta 

((lj)J B , Jm\T Jm \0) = ^{(lj)J B , J\\Tj\\0) (76) 
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and reducing the products of two rotation matrices and two spherical harmonics to linear 
combinations of spherical harmonics 



1 1 



vy*>(g)vw\sy = V^J2(- l ) 1/2+u+J fj /2) y 2 V I m) Yjm{§) (77) 

fa. ~V^r[-M M' M')[o J^'(P). 

(78) 

with fj^ = the Fano tensor for spin-1/2 polarization, we obtain 

BV'r = EEE E E E i l '-\-i)^ +v+J f { } ) ( j I j m \yjm{s) 

M b vv 1 JM lMjm p l'M'j'm' p J'M' V / 



x (-1)' * m »[f] y I, l M , _ J m , j (-l) Js J m yi> p m b - m > ) 

x (_i)i-h-™>U] ( 1 Z J ) (-i)J B -j-m ( J Jb J \ 
K J [Jl \ v M -m p J \m p M B -m J 

x((0')j B ,j'||r;i|o)*((Zi)j B ,j||r J ||o), (79) 

where we have transformed the Clebsch- Jordan to three-J coefficients. Next we perform 
the sums over third components of angular momenta in the above expression. Note that 
the total phase inside the sum can be simplified to 

phase = (-l)3+^(_i)^+i+''(_i)J-J'. (80) 
Therefore the following coefficient appears 



M B vv' MM 1 m p m' \ / 



I I' J' 

-M M' M' 



x 



\ V f ) f Jb J' ) \ I j \ j Jb J \ 
v' M -m' p J \m' p M B -m' J \ v M -m p ) \m p M B -m J 

(81) 
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We first perform the sum over v, z/, M, M' by using a 9-j coefficient 



1 1 



1/1/ MM 



* J 

' M 



I V J' 



3 



-M M' M' \ v M -rru 



\ V f 
v' M' -mi 



LM 



J J' L 
M M M 



L J 



M m p —m' p 



J J' L ) 

\ 1 J 
\ V f J 



> . 



(82) 



Next we compute the sum over m p , m' Mb using a 6-j coefficient 



m p m' p MB 



L 3 3 
M m p —m' p 



L J' J 

M —m! m 



3 Jb J 
m p M B —m 

r l f j] 

' 1 Jb 3 3' j 



J' 



m' p M B —wl 



(83) 



Then the S'-coefficient (81) results 



s = EW 2 (-i) 

LM 

' J J' L 

x < \ I j 

\ 2 1 3 



J+l'+j+J B +m' 



J' J 



L J' J 

Jb 3 f 



J J' L 



M -m! m ] \ M M' M 



(84) 



To finish we insert the result (84) into (79), and define indices a, a' corresponding to the 
quantum numbers of the final states 



a=(l,j,J), a' = (l',f,J') 



(85) 



and a coupling coefficient 



$ a , a (j,j',L) = v^w[ni>ii>-'][^][^][Ji[i](-i) ,+i+Ji,+L+j+j ' 

J J' L ) 



x 



I V J' 




L J' J 

Jb 3 f 



I 3 
I' 3' J 



(86) 



The final expression for B is 



B .J'J 



^EE E i l '-\-ir( 

Ij JJ'LM \ 



J J' L 
m —m! M 



<t> a , a (J,J',L) 

\ iii — iii \ i i 

[j [/j' J Jl 

x [Yj{s)YAv')],-m ((I'JVb, J'WTj'WOnmJs, J\\Tj\\0) , 



(87) 
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where the coupling between two spherical harmonics has been used 



[Yj(s)Yj,(p')] LM = 



J J' L 
M M' -M 



) 



Y JM (s)Yj, M ,(p') . 



(88) 



MM' 



Although the result given in eq. (87) is formally identical, with the exception of the 
factor 1/2, to the one obtained in [49] for the case of polarized nuclei, there exists a basic 
difference concerning the polarization coefficient ^ a i a (J, J', L), which contains all the 
information on the polarization properties of the particles in the initial and/or final state. 
Note that in the present case (spin-1/2 polarized particles), the angular momentum in the 
expansion of the rotation matrices, J7, only takes the values 0,1. The case J = is the 
only one surviving when the final nucleon is not polarized, i.e., when summing the cross 
sections for ±s values. In this case the present formalism reduces simply to the standard 
unpolarized one of ref. [41]. In fact, for J = we have J 1 = L and the reader can prove 
after some Racah algebra, that $ a ' a (0, L, L) reduces to the expression given in eq. (All) 
of Ref. [41] for the unpolarized case. 

Moreover, using the properties of the 9-j symbol, the following important symmetry 
property is found for the polarization coefficient under exchange of the indices 



This property coincides with the one already presented in [49] in the case of polarized 
targets. Since the multipole expansion of response functions performed in [49] was based 
on this symmetry, then the same formalism can be applied to the present case. In this 
way one arrives to eqs. (37-45) (see [49] for more details on the expansion). 
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